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Abstract 

This brief article treats question of fundamentality of the translates of a polyharmonic 
spline kernel (also known as a surface spline) in the space of continuous functions on a 
compact set C R'' when the translates are restricted to f2. Fundamentality is not hard 
to demonstrate when f2 = R'' or when fl is compact but a low degree polynomial may be 
added; the challenge of this problem stems from the presence of the boundary, for which 
all successful approximation schemes require an added polynomial. 

When f2 is the unit ball, we demonstrate that translates of polyharmonic splines are 
fundamental, by considering two related problems: the fundamentality in the space of 
functions vanishing at the boundary and fundamentality of the restricted kernel in the 
space of continuous function on the sphere. This gives rise to a new approximation scheme 
composed of two parts: one which approximates purely on dil, and a second part involving 
a shift invariant approximant of a function vanishing outside of a neighborhood n. 

1 Introduction 

Interpolation using the thin-plate spline radial basic function (/>(r) = r^logr requires the 
addition of a linear polynomial. This requirement is necessary to guarantee the solvability of 
the interpolation problem (see e.g. Light [7j). 

It can also be shown, that for any compact set Q C R'^ (this will always denote a compact 
set in this paper) functions of the form 

J2 ay(l){\x - y\)+Pi{x), 
yev 

where y C is a finite, are dense in C(Q). Here pi is a polynomial of degree 1. The uniform 
norm on is denoted by || ||oo,n- If the domain is the whole of Euclidean space we shall drop 
the n. 
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Definition 1. A set A is fundamental in normed vector space X if, for every x £ X , given 
e > 0, there exists a finite subset Y C A and coefficients Oy, y £Y , such that 



^ayy 

y& 



< e. 



In Cheney and Light [4j there is the excehent Chapter 18 which discusses the fundamen- 
tahty of translates of positive definite and conditionahy positive definite functions of order 
1. The arguments presented therein fohow that of Brown [3j, which are measure-theoretic 
and non-constructive. There is a measure theoretic proof for density of the thin-plate spline 



(j){r) 



log r in R , shown to us by Johnson [6j , but this approach does not generahse to 



other space dimensions. We do not follow that route here, but develop constructive proofs of 
fundamentality. 

In this paper we will show that translates of the polyharmonic splines 



r^^-'^, d odd, 
„2A:-djQg^^ d even. 



are fundamental in the continuous functions on the closed unit ball 23'^ = {x £ R*^ 



X 



<!}■ 



We shall denote these functions by C(23 ). As a step on the way we shall prove that for any 
compact set Q, C R"^, such functions are fundamental in the set Co{Q), of continuous functions 
which are zero on the boundary dQ. 



2 Density in Co{Q), C{dn) and C{Q) 

In this section we discuss approximation by integer shifts of a function satisfying a simplified 
type of Strang-Fix condition. We use this scheme to approximate functions vanishing on the 
boundary of a compact set Q. We then discuss that the polyharmonic B-splines - particularly 
that they satisfy such a condition. Finally, we demonstrate that the polyharmonic splines are 
fundamental in C(il) if dil, is such that their restriction is fundamental in C{dil,). 



2.1 Approximation in Euclidean space 

A standard way of constructing an approximation to a function is via convolution with an 
approximation to the identity. Given / S C(R'^) and g G Li(R'^), the convolution is 



f *g{x) = 




f{x - y)g{y)dy. 



If we further insist that 5 = 1; then we can construct the approximate identity gh{x) = 
h~'^g{h~^x), X G R"^, h> 0. In this case, gh* f ^ f uniformly as /i — > / G Co{Q). 

If we discretise the convolution f *gh on the /i-scaled multi-integer grid (this is sometimes 
referred to as the semi- discrete convolution) we arrive at the quasi-interpolant 

Sh{f,x) = ^ f{hz)g{h'^x - z). 
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We apply Sh to functions that are uniformly continuous on W^, i.e. functions for which 
the modulus of continuity uj{f, 6, := sup{/(x) — f{x — y) : x £ ^l, \\y\\ < 5} (suppressing Cl 
when Q, = R"^) satisfies 

co{f,6)=u{f,6,K'')^0, 
as (5 — )• 0. To prove this convergence we insist that g satisfies the following three conditions. 
Assumption 2. We assume the function g G C(M'^) satisfies the following: 

(I) The translates of g form a partition of unity: i.e., for any x G R'^, S^gZ'* di-^ — z) = 1. 
(II) There exists an M such that for every x G R'^, S^eZ'* ~ ^)\ ^ 
(III) There exists an M' such that for every x G R'^, X^^gz*^ ~ ^Ibl^^ ~ z)\ < ■ 

The rates of convergence of the quasi-interpolant operator Sh and other shift-invariant 
approximation schemes have been studied extensively under the umbrella of shift invariant 
approximation. The approximation power of such schemes are related to Strang - Fix condi- 
tions (conditions on behavior of the Fourier transform of g at 27rZ). We note that Assumption 
[2] guarantees that g satisfies a Strang - Fix condition of order zero (i.e., g vanishes at the 
27rZ \ {0} with a simple zero) since the semi-discrete convolution reproduces constants. See 
[5] for a discussion of approximation rates for this operator. We do not require the full force 
of this theory and develop what is necessary for our purposes in the following proposition. 

Proposition 3. Let f be hounded and uniformly continuous on R'^. Then, given e > 0, there 
exists an H > such that \\f — Sh{f, ■)\\od < e for all h < H . 

Proof: Since the translates of g form a partition of unity (item (I) of Assumption [2]) , 



\f{x)-Shif,x)\ 



X — z, 



zGZd 



< 5^ |/(x)-/(M||5(/i-^ 



X — z) 



X — z 



< E \fi^)-fihz)Mh-' 

hzGBg (x) 

+ E \fix)-fihz)\\g{h-'x-z 

hz^Bs{x) 

< u;U\6) E \gih-'x-z)\+2 

hzGBg (x) 



oo E 1^*^^ ^X-Z)\, 
hz^Bs{x) 



where Bs{x) = {?/ G R'^ '■ \x — y\ < 6} is the ball of radius 6 centred at x. Let us choose 6 
such that uj{f,6) < e/(2M), and following this choice of 5, choose h < e5/(4M'||/||oo)- Then, 
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using Assumption [2] (II) and (III), 



\fix)-s,,{f,x)\ 

< Ma;(/,(5) + 2||/||oo sup \h-^x-z\ \h-^x - z\\g{h''^x - z)\ 

< l + 2||/|U^M'<| + ^ = e. ■ 

2.2 Approximating functions vanishing on the boundary of a compact set 

The difficulties one faces when proving density results via discretisation of convolutions happen 
due to the boundaries of the region. Because of this it is much easier to prove results where 
the function is zero on the boundary of the domain of interest. This shall be the case for this 
section. Since our target function is continuous on the compact domain Q, it is uniformly 
continuous. If we let / to be defined on the whole of R"^ by setting its value to zero on the 
complement of Q, it is easy to show that / is also uniformly continuous on R'^. 
We define our approximation 

s^if,x):= Yl f{hz)g{h-^x-z), xeR-^. 

Our first result, which we will use later, says that if / is small on O then s^(/, •) is small 
everywhere. 

Lemma 4. Suppose f G C(il) and ||/||oo,r2 < £■ Suppose further that g satisfies item (II) of 
AssumptionlM Then, ||s^(/, OIloo < Re- 
proof: Since g satisfies Assumption [2] item (II) we have, for x G R'^, 



Y f{hz)g{h-^x - z) 
hzen 



< ||/||oo,n Yl \9ih-'x-z)\<Me. 
hzen 



Since, for x ^ Q, s^{f,x) = Sh{f,x) as a corollary of Proposition [3l we get our first 
approximation result for a compact domain: 

Corollary 5. let f G Cq{0,) and g satisfy Assumptions^ Then, given e > there exists an 
H > such that \\f - s^(/, ■)\\oo,n < e for all h < H. 



2.3 Poly harmonic B-sphnes 

We employ the fact that, in the distributional sense, the iterated Laplacian A^cjid^k = Cd^k^ 
for some real number Cd^k- Let {ej : j = 1, • • • , d} be an orthonormal basis for R'^. We can 
now form the discrete Laplacian of a continuous function / 

d 

2df{x)-Y{fi^ + hej) + f{x-hej)) . 



A/^/(x) 



1 

h^ 
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Then, as suggested by Rabut ^0], we can form the polyharmonic B-spline^ 

1 



B. 



d,k,h 



-A 



When h = 1, we drop the third subscript: B^f^ := B^f^^i. It follows that B^i^ = 
XlzeZ'* ^z4'{\x — z\) for a sequence {cz)z^z''- satisfying = when Yl'j=i > ^- rescaling 
the coefficients, Bd,k,h = Y.zeZ'i '^z4>{\x -hz\). 



Proposition 6 (From [lOJ). The function B^ k = -Brf,fc,i satisfies Assumption\^ Additionally, 
the functions satisfy the dilation condition: for any h > 0, 



{Bd,k)h {^) = h-'^Bd,! 



Bd,k,h{x). 



Proof: Assumption [2] follows from [10]. Specifically, / follows from [101 Theorem 4] 
while (II) and (III) are a consequence of the decay \B(i}^{x)\ = 0{\x\~'^~'^) in [lOJ Theorem 2 

(i)]- 



When d is even, (l){\hy\) = h 



2k~ 



y\) + {h"'-'' log h)\y 



2k-d 



SO 



E c.//.<^fb - z\) = h-^^ 



zehZ'> 



hx — /i^ 



since X^^eZ'* 1^ ~ ^|2fc-d _ q_ \Y]^gj^ ^ jg q^j^^ (^^^^ homogeneous, and the dilation condition 
follows by a similar (in this case simpler) argument. ■ 

Although we are primarily interested in proving density properties of -Bd.fc? there exists an 
extensive literature on convergence rates for the operator Sh and related schemes. We note in 
particular that the function i?^ ^ satisfies a Strang-Fix condition of order 2m, which can be 
observed directly from the Fourier transform of B^^k in [10, Theorem 2]. We refer the reader 
to [5], Theorem 3.6] for a discussion of the approximation orders of this kind of approximation 
scheme. 



2.4 Approximating on compact sets 

We now give our main result. 

Theorem 7. Let 0, be a compact set. Assume that the restriction {x,a) i— )• (f)d,ki\x — a\) is 
fundamental in C{d^}). Let f G C{0,). Given e > there is a finite subset Y C 0, and a set 
of coefficients {oy : y £ Y} such that 

y&Y 

for all X £ 0,. 



^ Rabut 's definition differs from tliis one by a factor of h'^. 
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Proof: In the proof of this theorem we will assume that the kernel g of Assumption [2] 
is Bd^k- By assumption, we can select a finite subset W C dQ and coefficients {uy : y £ W} 
such that 



< 



4M' 



fix) - J2 '^y^d,k{\x - y\ 
yew 

for all X G dQ. Here M is the specific constant used in Assumption [2] item (II). 
For 7 > let us define the following neighborhoods: 



a 7 neighborhood of 0: 



Cl.y := {x I dist(x,r2) < 7}, 



• an annular neighborhood of the boundary: 

:= {x I dist(x,c?r2) < 7}, 

• a one-sided neighborhood 

S^:={x £n \ dist(x, dn) < 7}. 

Now let us set r{x) = f{x) — J2y£W ^y'^<i,k{\x — y\) , x £ Q. Then, since r is continuous and 
small on 50, there is a (5 > such that 



\r(x)\ < 



2M 

for all x G A^. By multiplying by a continuous cut-off function a 



a{x) = < 



dist(a;,9r2) 







■g , < dist(a;,r2) < 5 

dist(x, Q) > 5 



there is a continuous function r = a x r such that r = r on J7, r(x) = when x fi^, and 
lkl|oo,Sj < ll'^lloojAa < e/(2M). Thus, using Lemma|l]we have, for any h > 0, 



stir) 



e 

< -. 

- 2 



Since B^k satisfies Assumption [2] there exists an H > such that, for h < H 



r — s^^{r) 



< -. 

00 2 



(1) 



(2) 



We approximate f{x) with s{x) := "^y^w ay(j)d,k{\x — y\) +8^"^^^ {r,x). Then, since 
Sh'ir) - st'{r), 

||/-s||oo,Q = \\f - ^ ay(pd,k{\x - y\) - s^'^^^{r)\\oo,n 
yew 

< \\r-s]^^{f) + st'{r)\\oo<e, 



n\Ss, 
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for h< H, using (P and ([2]). 

Finally, by Proposition El the terms making up the approximant s^'^^^ {f) can be expressed 
in terms of the hZi'^ shifts of (f). Indeed, we have 

h~'^Bd,k{h~^x- z) = Bd,k,h{x-hz) = h-^^^CQ(t){\x-hz-hC\) 

with = when X]j=i Kjl > ^- This means that so long as h < 6/k the centers used to 
construct s lie in 17. To whit, consider the hz term in s^'^^^ (namely the term B(i^k{h~^x — z) 
centered at hz). It uses centers inside the closed ball B{hz, kh). If there is y € B{hz, kh) so 
that y ^ i}, then there is a point on the boundary between y and hz (and, hence, hz is nearer 
to the boundary than 5) . ■ 



3 Density in C(:B^) 

In this section, we apply the scheme developed in Section 2 to prove fundamentality in 17 = 'B'^. 
We will first show that we can approximate arbitrary continuous functions on the sphere using 
linear combinations of polyharmonic splines centred on the sphere. This fact, in conjunction 
with Theorem [7] will deliver our density result. 

3.1 Approximation on the unit sphere 

We begin by introducing the Fourier- Gegenbauer expansion of a function -0 : [—1,1] — )• R- 
The Gegenbauer polynomials P^, A > —1/2, k = 0,1,..., are orthonormal polynomials (we 
will assume this normalisation in what follows) on the interval (—1,1) with respect to the 
weight w^{t) = (1 — t'^)^~'^^'^. In other words, 

j ^ P^{t)Pj^{t)w\t)dt = 5kj, k,j = 0,l,.... 

The Fourier- Gegenbauer coefficients of the function / : [—1, 1] — )■ R are 

«fc(/) = fjPk\t)w\t)dt. (3) 

The first step of our proof is to approximate / G C(23'^) on the boundary S'^"^ using 
polyharmonic splines, li x,y G S''"^ then |x| = \y\ = 1 and we can rewrite 4'd,k{\x — y\) = 
'>Pd,kix ■ y), where 

, . _ J (2 - 2r)^-^/2^ d odd, 

^d,fc(r-j - I 1 ^2 _ 2r)'=-'^/2 log(2 - 2r), d even. 

A crucial result for us is found in Sun and Cheney }111 Theorem 2]. 
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Proposition 8. Fix d and let X = {d — 2)/2. Let -0 G C[— 1, 1]. In order that the the set 
{x I— 7- iIj{x ■ y) : y G S'^"^} be fundamental in C(S'^~"^) it is necessary and sufficient that 
a^{f)^0, k = 0,l,...- 

In j2] the Fourier-Gegenbauer coefficients for ■0^^^ are computed. For our immediate goal, 
it is not important what the values are, only that they are nonzero (see below for more details 
about a^(/)). We state this without proof here. 

Lemma 9. For d > 1 and k > d/2, The Fourier-Gegenbauer coefficients of (l)d,k cif^ o-ll 
nonzero. 

An immediate consequence of the above lemma, using Proposition — [Sj we get the first new 
result presented in this paper: 

Theorem 10. The set {(j)d,ki\ " -^1) • V ^ ^'^'^} dense in C(S'^~^). 
This, in conjunction with Theorem [7] proves the following. 

Theorem 11. Let f £ C(!B'^). Given e > there is a finite subset Y C 23"^ and a set of 
coefficients {ay : y G 1"} such that 

\fix) - ^ ay(l)d,ki\^ - y\)\ < e 

for all xe'B'^. 



3.2 Approximation on balls of arbitrary radius 

We now consider density of (j)d,k in ^"^(p)) the ball of radius p > 0. In this case, the proof of 
Theorem [7] applies, if the analog to Theorem 1101 is established for the sphere of radius p > 0: 
§<i-^{p) ■={xeR'^\\x\= p}. 

Equivalently, this problem can be considered by simply rescaling ct>d,k{px) = p'^'^^'^4>d,k{x) 
and considering the density of (pd^k on 'B'^ ~ clearly, the translates of (/)d^k dense in the set 
'B'^ip) if and only if the translates of (pd^k are dense in 'B'^. A quick computation shows that 

T JkP*'"'^, d is odd 

\\x\ log |x| + log/>|x| , a IS even. 

In odd dimensions, (pd,k = (t>d,k is scale invariant, while in even dimensions it is necessary to 
add a polynomial term (in this case, 4)d^k depends on p as well as d and k). Thus, the odd 

dimensional case follows directly from Theorem [11] and we consider d G 2N. The operator 


A annihilates polynomials of degree 2/c — 1, so the only hinderance to obtaining density in 
S'^(p), is the density of 0d^fc|gd-ixs<i-i restricted to S"^^^, (and it is only necessary to consider 
this when d is even). 
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The restriction of (f)d,k to the unit sphere is 

4>d,k{^-y) = \{2-2x-yf-'''^\og{2-2x-y) 
+ \og{p){2-2x.yf-'''^ 

Since d is even, p{t) = {2—2t)^~'^/'^ = X]j=o^^ ^/(*) ^ polynomial. Its Fourier-Gegenbauer 
coefficients are computed in (2.12)]. For A = they are 

(-1)^-22^ ( d \ fd \ ik-d/2)\T{k-l/2) 
' 40F V 2 y V2 J-{k-d/2-j)\{k+j + d/2-2y. 

Fourier-Gegenbauer coefficients of V'd.fc are computed in [2, (2.19)]. For j < k — d/2, we have 
that the coefficients (with A = ^) of ipd,k{t) = Ejlo ^/(*) satisfy 6^ = x rfc^^^j with 

tmj = ^'(A:-d/2 + l) + ^(A;-l/2) + 21og2 

-^{k -d/2 + 1- j) - ^(A; + j + d/2 - 1) 

k-d/2 k-1 k-d/2-j k+d/2+j 

EL ^ — > Z ^ — > i ^ — > 1 

n ^ 2n-\ ^ n ^ n 

n=l ra=l n=l n=l 

Here ^' is the digamma function F'/F and ^(m + 1) = -7 + Y!^,=\ ^ (cf- [T, (6.3.2)]) while 
^(m + l/2) = -7-21og2 + E"=i^ (cf. [1,(6.3.4)]). 

It follows that when p > is chosen so that log/) = —Tk^j for some j < k — d/2, then the 
restriction of the kernel (pa^k to S^^^(p) is not dense. 

Example 1. We consider the thin plate spline example k = 2,d = 2. In this case, degenerate 
examples occur when log/? = — T2,2,i or — T2,2,o- ^ simple computation shows that 

1 1, 12-6-3-2 1 
r2,2,o = l + 2-l-(l + - + -) = = - 

and 

111, 1 
r2,2,i = l + 2-(l + - + - + -) = -- 

from which we obtain p = e^l^ and p = e~^/^'^. 



4 Conclusion 

The results of this paper are a first attempt at treating a challenging problem in RBF theory: 

Identify radial functions </> so that {(/){■ - C) | ^ G 23^^} is fundamental in C(23'^). (4) 

This can be viewed as a modification of the density problem considered by Pinkus [9j to 
characterize the radial functions whose translates (from R*^) form a fundamental set in C{K) 
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for a compact K. It has been shown that for any radial 4> G Li(R'^) H C(R'^), the set of 
translates {(/>(• — ^) | ^ G I^'^)} is dense in C{K) for any compact K. This is not exhaustive: 
the class of functions for which this holds is significantly larger. 

The solution of problem may be a much smaller class of functions, simply because 
it modifies the original problem in a nontrivial way. Making a challenging mathematical 
problem more difficult does not generally make it more elegant or interesting. We do so 
in spite of this observation because the resulting problem has value also in understanding 
of the role of the boundary in RBF approximation. It is well known that the presence of a 
boundary - where centers are kept to one side - is detrimental in a number of ways (it curtails 
approximation rates, forces instability by way of increased Lebesgue constants, causes blow 
up of error near the boundary, etc.). We pose this question because we wish to know when 
(or if) a polynomial term is needed to construct a good approximation and, secondly, in the 
hope that it will generate approximation schemes that adequately treat the boundary effects. 

As a final remark, we consider the approach to solving ^ suggested in this paper in more 
general light. Of the two components of the approximant - the part that approximates in 
Co(ri) and the part that approximates in C(50) - we have little to say about the former, since 
few radial basis functions are known to have a dilation invariant, rapidly decaying localization 
(this is a unique perk of working with polyharmonic splines). The latter component can be 
generalized quite easily, however. 

For a radial function G C(M'^), the restriction of (x, y) i— )• (f){x — y) to S"^"-*^ satisfies the 
relationship (j){x — y) = "0(1 — x ■ y) for some univariate function -0 : [0, 2] — t- M. The Fourier- 
Gegenbauer coefficients of 1 1— ?■ 'tp{l — t) are tied to the (distributional) Fourier transform of (j). 
In particular, if (p G -Li(M'^) has a Fourier transform in Li, then aj" = f^-^ ip{l — t)Pj'{t)w^{t)dt 
is determined by the Weber-Schafheitlin type integral: 

a]= / ta{t){Jj+x{t)fdt (5) 
Jo 

for all j (see Proposition 1] for a proof). Here A = and a : [0,oo) — )■ M satisfies 
cr(l^l) = 4>{£,) and is a first order Bessel function. It follows that if io"(t) ( Jj+A(i))^df is 
nonzero for all j (for instance, if ait) > and is nonzero on a set of positive measure) then 
the restriction of </> is dense in S*^"^. 

Recently, zu Castell and Filbir fT2j and independently Narcowich, Sun and Ward [8] have 
extended this idea to include a class of radial functions that exhibit polynomial growth. In 
|12l (4.5)] and in [U Proposition 3.1] a similar relationship to ([5]) is shown for conditionally 
positive definite (CPD) function of order k: Fourier-Gegenbauer coefficients 

poo 

Jo 

for j > 2k + 1. Here the distributional Fourier transform / restricted to M'^ \ {0} (one may 
use the so-called generalized Fourier transform which is necessarily radially symmetric), is 
represented by the positive Borel measure a: i.e., = /o°° Xjes^-i V'(^'^) do; d(T(r) for 

Schwartz functions ip vanishing at the origin like 0(|'^P'^). 
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In this case, if for instance a is represented by a function that is positive on a set of positive 
Lebesgue measure, then modulo addition of a polynomial of degree 2/c + 1, the restriction of 
/ is fundamental in S''"^. This is more or less what happens in Section 3.2, where the 
two condtionally positive definite functions ^^.k f'^id (j)^^}; have the same generalized Fourier 
transform ^ ICI"^*^) but their Fourier-Gegenbauer expansions disagree at low indices. 
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